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ABSTRACT 

We present an analytical formalism for the description of Blazhko RR Lyrae light curves in 
which employ a treatment for the amplitude and frequency modulations in a manner similar to 
the theory of electronic signal transmitting. We assume monoperiodic RR Lyrae light curves 
as carrier waves and modulate their amplitude (AM), frequency (FM), phase (PM), and as 
a general case we discuss simultaneous AM and FM. The main advantages of this handling 
are the following: (i) The mathematical formalism naturally explains numerous light curve 
characteristics found in Blazhko RR Lyrae stars such as mean brightness variations, compli- 
cated envelope curves, non-sinusoidal frequency variations, (ii) Our description also explains 
properties of the Fourier spectra such as apparent higher-order multiplets, amplitude distribu- 
tion of the side peaks, the appearance of the modulation frequency itself and its harmonics. In 
addition, comparing to the traditional method, our light curve solutions reduce the number of 
necessary parameters. This formalism can be applied to any type of modulated light curves, 
not just for Blazhko RR Lyrae stars. 

Key words: methods: analytical — methods: data analysis — stars: oscillations — stars: 
variables: general — stars: variables: RR Lyrae 



1 INTRODUCTION 



The Blazhko effect l |Blazhko"||1907^ is a periodic amplitude and 
phase variation of the RR Lyrae variable stars' light curve. The typ- 
ical cycle lengths of these variations are about of 10-100 times 
longer than the main pulsation periods (0.3—0.7 d). Almost half 
of the RR Lyrae stars pulsating in their fundamental mode (type 
RRab) and a smaller but non-negligible fraction of the first over- 
tone mode pulsating stars (type RRc) show the effect jjurcsik et 



al. |2009c[ IChadid et aL 2009 , Kol enberg et al. |2010[ |Benko et 



al. ||2010' i. It is usually interpreted as a modulation or a beating 



phenomenon, but both hypotheses have their own problems. The 
beating picture describes the main feature of the light curves and 
Fourier spectra well (see Breger & Kolenberg 2006| [Kolenberg et| 
|al. ||2006l >, but reproducing phase variations, multiplet structures 
found in certain stars' Fourier spectra jlurcsik et al. |200¥{ IChadid] 
|et al. |2010| l in this framework is not possible. On the other side, 
the stars showing doublet stmctures in their Fourier spectra (|AF] 
[cock et al. |200 2003, Moskalik & Poretti |2003,) seemed to be 
contradicted with the modulation picture. 

In this paper we describe the Blazhko effect as a modulation, 
and derive the mathematical consequences of this assumption by 
developing a consistent analytical framework. Using this frame- 
work we demonstrate that many light curve characteristics are nat- 
urally identified as mathematical consequences of the modulation 
assumption. By disentangling these features we get closer to the 
physics behind the Blazhko effect. 
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The possibility of the modulation/Blazhko effect have been 
raised for many types of pulsating stars from Cepheids to 5 Scuti 
stars (see e.g. Koen [2 001 , Henry, Fekel & Henry 2005; Moska^ 
|lik & Kotaczkowski [2009,.Breger .20 10, .Poretti et al. .2011,) . The 
main motivation of this paper is the mathematical description of the 
Blazhko RR Lyrae stars' light curves and investigate their proper- 
ties. Most of our results can be applied directly to any other types 
of variable stars, where modulation is suspected. Our deduced for- 
mulae and the related phenomena may help to prove or reject the 
modulation hypothesis. 

The basic idea of this paper was raised in'Benko et al. mQQ9^. 
Modulation is a technique that has been used in electronic commu- 
nication for a long time, mostly for transmitting information signal 
via a radio carrier wave. In those cases, the carrier wave is a sinu- 
soidal electromagnetic (radio) wave that is modulated by a (gener- 
ally non-periodic) information signal (e.g. speech, music). In this 
paper the formalism developed by engineers for broadcasting radio 
signals has been modified such a way that we assumed a monope- 
riodic non-modulated RR Lyrae light variation as a carrier wave. 
While for communication usually only one type of modulation is 
applied we allow both types of modulations (amplitude and angle). 

In Section |2] we present a collection of classical formulae 
that are well-known in physics of telecommunications ( [Carson | 
[19221 [van der Pol [[l930l [Roder [[T93l] >. Some of the more com- 
plicated cases (multiple modulations, recursive or cascade modula- 
tions) were investigated by mathematicians who developed the the- 
ory of electric sound synthesizers in the years of 1960s and '70s. 
The formulae are modified to describe the modulated light curves in 
Sec.|3]and investigated by a step-by-step process from the simplest 
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cases to the more complex ones. Section 4 compares the numerical 
behaviour of the traditional method and this one. Section 5 summa- 
rizes our results. 



2 BASIC FORMULAE 

In this section we briefly review some classical definitions and for- 
mulae (see e.g. |Newkirk & Karlquist |2004[|Schottstaedt |2003l l that 
will be used through the next sections. The simplest periodic signal 
is a sinusoidal function. It has three basic parameters: amplitude, 
frequency and phase and any of these can be modulated. 



2.1 Amplitude modulation 

The amplitude modulation (AM) is the simplest of the three cases. 
Let the carrier wave c{t) be a simple sinusoidal signal as 



c(t) = Uc sin(27r/ct + (pc), 



(1) 



where Uc, fc and ipc constant parameters are the amplitude, fre- 
quency and initial phase of the caiTier wave, respectively. 

Let Uni{t) represent an arbitrary waveform that is the mes- 
sage to be transmitted. The transmitter uses the information signal 
Uin{t) to vary the amplitude of the carrier Uc to produce a modu- 
lated signal: 

(7am (t) = [Uc + Um (t)] sin(27r/cf + ipc). (2) 
In the simplest case, where the modulation is also sinusoidal 

Um{t) = U^ sin(27r/„,t + if^,). (3) 

Substituting Eq. Jsjl into |2| and using basic trigonometrical identi- 
ties, expression O can be rewritten as 

UAM{t) = Uc Sm{2Tvfct + ifc) + 

r/A 

{sin [27r (/c - t + p-]+ sin [2tt {fc + fm) t + ip+]} , 



(4) 



where(y9^ = ipc — (p^ + iT/2 and (p'^ = i/^c + V^m ^■''"/2. The initial 
shifts (±7r/2) appear because we used (through all this paper) a 
sinusoidal representation instead of sin and cos functions. 

The exact analytical Fourier transformation of (|4| is given in 
Appendix [A[ however, the basic structure of the frequency spec- 
trum can be easily read off from Eq. Since the Fourier spectrum 
of a single sinusoidal function shows a peak at the frequency of 
the sinusoid, from the above expression |4]l the well-known triplet 
structure composed of the peaks fc and fc ± /m can be seen. The 
amplitude of the side-peaks fc ± fm are always equal. The Fourier 
amplitude of the carrier wave (7r\/27r[/c), that represents the energy 
at the carrier frequency, is constant. 

The ratio of the carrier wave amplitude A{fc) and the side 
peaks A{fc ± /m) are connected to the modulation depth. We 
rewrite Eq. ([2]l as 



(7am (t) 



1 + 



Un. (t) 



Uc 



(5) 



If Uin(t) is a bounded function, let U^^^ represent the maximum 
value of this modulation function, then modulation depth is de- 
fined as h = Um^^ /Uc- In the above discussed sinusoidal case 
h = Um/Uc and A{fc ± /m)M(/c) = In other words, the 
amplitude of the central peak is twice of the side peak highs. 



2.2 Angle modulations 

The phase and frequency modulations together are called angle 
modulations. Since, when we assume the sinusoidal carrier wave 
Eq. ^ as c{t) = Uc sin[9(t)], the Q{t) = 2Tvfct + ipc denotes the 
angle part of the function. 

Phase modulation (PM) changes the phase angle of the carrier 
signal. Suppose that the modulating or message signal is (7m (t), 
then e{t) = 2-Kfct + [ipc + Um{t)]. Let Uui{t) be again a 
bounded function. In this case we can define a constant as: fcpM ~ 
\Um'^^{t)\/2. This transforms the modulated signal 



U-puit) = Uc sin ^2Tvfct + k-puUm (t) + <^c] 



(6) 



where \ Um{t)\ ^ 1. The instantaneous frequency of the modulated 
signal is 



fit) 



de 



fc + fcpi 



(7) 



Frequency modulation (FM) uses the modulation signal 
Uin{t) to vary the carrier frequency. Q(t) — 2nf{t)t + ipc and 
here the instantaneous frequency f{t) is modulated by the signal 
of fcpM(7m(t) as 



fit)=fc + kFMUUt)- 



(8) 



In this equation fcpM is the frequency deviation, which represents 
the maximum shift from fc in one direction, assuming U^it) is 
limited to the range (—1, . . . , +1). Using the definitions of the in- 
stantaneous frequency and phase, expression ^ can be rewritten as 

Q{t) = 2-Kfct + 27rfcFM /q* Um{T)dT + ific- The modulated signal 



2Tvfct + 2nkFM / Um (t) dr + ipc 



(9) 



Ufm (t) = Uc sin 



This definition of FM is the least intuitive of the three Eqs |2] |6] 
and|9j. If we compare Eqs |6| and |9]( we realize that the modu- 
lation signals are in derivative-integral connection with each other. 
In practice, we have modulating signals that can be represented by 
analytical functions. Therefore, when we detect an FM or PM sig- 
nal without any previous knowledge about them, it is impossible to 
distinguish between FM and PM signals. 

First, let the modulating signal be represented by a sinusoidal 
wave with a frequency fm - The integral of such a signal is 



u" 

IT 



sin ( 2-Kfmt + ip„ 



(10) 



fc] , (11) 



27r/n: 

Thus, in this case Eq. l|9j gives 

UFM{t) = Uc sin |^27r/ct + J?sin (2-Kfmt + < 

where the modulation index is defined as r; = (kp-MUm)/ fm- 
Eq. ([TTJ can be deduced from Eq. |6| as well. The only difference 
is the value of 77 = kpyiUm, which is independent of the modu- 
lation frequency fm- Let us transcribe Eq. ijTTJ using relations for 
trigonometrical and Bessel functions l |Abramowitz & Stegun |1972l > 



Ufm (t) = Uc ^ Jk iv) sin [27r {fc + kfm) t + kipm + 



^c\ 



(12) 



where Jk {rj) is the Bessel function of first kind with integer order 
k for the value of r/ (Fig.[T|(; ipm denotes either (pm or (pm- This 
formula is known as the Chowning relation l |Chowning [1973^ . Al- 
though it had been deduced by many different authors formerly. 
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Figure 1. The graph of the first three and the 10th Bessel functions of first 
kind with integer order k for the value of x. Jq {x) - (red) continuous line, 
J\{x) - (green) dashed line, J2{x) - (blue) dotted line, Jio(a:) - (purple) 
dash-dotted line. 



Chowning recognized its key role in the electronic sound creation 
method called FM synthesis. 

Similarly to Eq. |4| expression ^12) also helps to imagine the 
spectrum. It is made up of a carrier at fc and symmetrically placed 
side peaks separated by /m. The amplitudes follow the Bessel func- 
tions. The behavior of the Bessel functions is well known: except 
for small arguments (x < fc), they behave like damped sine func- 
tions (see also, Fig. [TJ. For higher indices the higher order side 
peaks gradually become more and more important. As a conse- 
quence, the amplitude of the central peak gets reduced. The fre- 
quency spectrum of an actual FM signal has an infinite number of 
side peak components, although they become negligibly small be- 
yond a point. 

If \rj\ < 1, we find that Jo(r;) 1, | J±i| = ri/2 and Jk ^ 
for fc > 1. That is, the spectrum can be approximated with an 
equidistant triplet similarly to AM, but the character of the sig- 
nal differs from AM: the total amplitude of the modulated wave 
remains constant. When ry increases the amplitude of side peaks 
also increases, but the Fourier amplitude of the carrier decreases. 
In other words, the side peaks could be larger than the central peak, 
on the other hand higher order side frequencies could also be of 
larger amplitude than the lower order ones. 

A more general case is formulated by |Schottstaedt |(l977[ l 



f/FM(i) = Uc sin 



+ J2 ' Sin (2n f^h + + 



kp — — oo ki — — oo 



1 



27r/ct + 



J2 kp {27V f^^h + <^i?') + <^c 



(13) 



Here the modulating signal is assumed to be a linear combination 
of a finite number of sinusoidal functions with arbitrary frequen- 
cies fm \ amplitudes (7nf' and phases i^m', p ~ 1,2, ... ,q. The 
spectrum contains equidistant frequencies on both sides of the car- 
rier frequency. The amplitudes of the side peaks fc ± kpf^^ are 
determined by products of Bessel functions. 



2.3 Combined modulation 

In practice, the electronic circuits that generate modulated signals 
generally produce a mixture of amplitude and angle modulations. 
This combined modulation is disliked in radio techniques but wel- 
comed in sound synthesis and as we will see they appear in the case 
of Blazhko RR Lyrae stars, as well. Let us overview the basic phe- 
nomena of combined modulations following ICartianiT] ( | 1 966) 1. We 
start with the simplest case: both the AM and FM are sinusoidal 
and their frequencies are the same. 

f7comb(i) = (7c (1 + ftsin27r/mt) sin [27r/ct+ 

r]sm{2'Kf^t + <l),-n) + ipc]- (14) 

By suitable choice of the starting epoch, without any restriction 
of the general validity we can set ipc ~ 0. Here 0^ is the rela- 
tive phase difference between the modulating FM and AM signals. 
Other designations are the same as before. The third term of the 
product l |14[ l is the same as in Eq. l |l 1[ ). therefore, after applying the 
Chowning relation l|T2j 



Ucomb{t) = Uc{l + ftsin27r/mt) • 

oo 

Jk iv) sin [27r {fc + fc/„0 t + , 



(15) 



This expression results in an infinite number of amplitude modu- 
lated waves. After trigonometrical transformations we get: 



(7comb (t) = Uc ('') si'i ifc + kfu,)t + kcf)^] + 

fc — — OO 

^ Jfe_i(?7) sin [27r {fc + kf„,) t + {k - 1) 4>,-n - |] + 
^ Jfc+i(j?) sin \2tt {fc + kf,^) t+{k+l) + f ] } • (16) 



It can be seen that each terms consists of three sinusoidal functions 
with different phases. On the basis of expression l |16^ , the spectrum 
of combined modulation ( |14^ is comprehensible as a combination 
of three FM spectra. The peaks are at the same places as the fre- 
quencies of the spectrum of l |12^ , but the amplitudes of a pair of 
side peaks are generally asymmetrical. Using some trigonometrical 
identities, the rules of summation of parallel harmonic oscillations 
and relations for Bessel functions we arrive to the expression for 
the Fourier amplitudes of a certain frequency: 



A{fc + kf^) ~ Uc { Ji{ft) ( 1 - ^ sin< 



— COS^(t>n,[Jk + l{v)-Jk-l{v)f> , (IV) 

(fc = 0, ±1, ±2, . . . ). I ntroducing the pow er difference of the side 



peaks as it was done by 



Szeidl & Jurcsik 



A'{fc + 



Ifm) — A {fc — Ifm), where I = 1,2,3 . . . and taking into account 
formula l |17| l, we get 



V 



(18) 



This formula is a direct generalisation of the formulae given by 
ISzeidl & Jurcsik | p00 9 1 for / = 1 and / = 2. It is evident, that 
this asynmietry parameter depends only on 0ni, the relative phase 
of AM and FM. The left hand side peaks are higher than the right 
hand side ones (A; < 0), if < (jim < t, otherwise the situation is 
opposite: n < < 2n and A; > 0. In those very special cases. 
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where i^m = or 0m ~ tt the side peaks' amphtudes are equal. 
We have to note that if one of the modulations from AM and FM 
dominates in the combined case (h <^ rj oi <^ 1), the amplitude 
of the side peaks are almost the same. 



3 BLAZHKO MODULATION 

RR Lyrae light curves traditionally are described by a Fourier series 
of a limited number of terms. In Blazhko modulated RR Lyrae stars 
the sum builds up from terms of harmonics of the main pulsation 
frequency, side peaks due to the modulation and the modulation 
frequency and even its harmonics: 



m{t) = Ao + ^A^ sin [2TvFit + 



(19) 




10 12 
Frequency [d"'] 



where either Fi = jfo, (j = 1,2, ...,n); or Fi — kfm, 
(k = l,2,...,m); or = j' fo + k' UAj' = l,2,...,n', 
fc' = 1, 2, . . . , m'); F, = j"fo - k"U, (j" = 1,2,..., n", k" = 
1,2,..., m"); /o and /m are the main pulsation frequency and the 
modulation one, respectively. The amplitudes Ai and phases "l>i 
are considered as independent quantities and determined by a non- 
linear fit. The necessary number of parameters for a complete light 
curve solution is 2N + 3 (amplitudes and phases and two frequen- 
cies and the zero point Ag). The number of parameters can be as 
high as 500-600 for a long time series of good quality (see e.g. 
IChadidetal. |2010p . 

In the next subsections we show how the modulation paradigm 
can be applied and what advantages it has compared to this tradi- 
tional handling Eq. \19\ . 



3.1 Blazhko stars with AM 

To start with, we discuss Blazhko stars' light curves with pure AM 
effect, although the recent space-born data suggest that all Blazhko 
RR Lyrae stars show amplitude modulation and simultaneous pe- 
riod changes ( |Chadid et al. |2010[[Benk6 et al. |2010[|Poretti et al. | 
|2010 l. We follow a step-by-step generalization process that allows 
us to separate effects more clearly. We note that the most striking 
feature of a Blazhko RR Lyrae light curve is the amplitude varia- 
tion, which is generally easy to find and in many cases the only 
detectable modulation (see |Stothers |2010l and references therein). 

To apply the framework described in Sec. |2.1| to an RR Lyrae 
light curve the course-book formulae need some extensions. 
We choose a continuous, infinite, periodic function with a non- 
modulated RR Lyrae shape as a "carrier wave". This function is 
described by the frequency fo and its harmonics, that is c*(t) : — 
m{t) if Fi = jfo in Eq. 
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Although, the exact analytical Fourier spectra of any of the 
modulated signals discussed in this paper can be calculated with- 
out any problems, at least in theory (see also Appendix|A|, to illus- 
trate the different formulae synthetic light curves and their Fourier 
spectra were also generated and plotted. An artificial light curve 
was constructed as a carrier wave with typical RR Lyrae parame- 
ters (fo = 2 d~^ and its 9 harmonics) on a 100 day long time span 
sampled by 5 min (insert in Fig.[2|. The Fourier transform of such 
a signal is well-known (Fig.[2](: it consists of the transformed si- 
nusoidal components given in Eq. |A2| (More precisely, due to the 
finite length of the data set and its sampling, the Fourier transfor- 
mations should always be multiplied by the Fourier transform of 
the appropriate window function.) 



Figure 2. Fourier amplitude spectrum of the artificial RR Lyrae light curve 
acting as a carrier wave of all modulated light curves constructed in this 
paper (main panel), and a part of the light curve itself (insert). 



Substituting c*(t) carrier wave into the definition of AM in 
Eq.(|5j: 



C/c* 



ao + ^ flj sin(27rj/ot + <^j) 



(20) 



Expression \20\ describes a general amplitude modulated RR Lyrae 
light curve, U*^^{t) is the modulation signal, U* is the amplitude 
of the non-modulated light curve. On the one hand, the non-zero 
constant term of ao is obligatory from mathematical point of view, 
otherwise the Fourier sum does not compose a complete set of 
functions. On the other hand, this value represents the difference 
between the magnitude and intensity means. More precisely, either 
we use physical quantities (viz. positive definite fluxes) or we trans- 
form normalized fluxes into magnitude scale. In this latter case the 
average of the transformed light curve differs from zero. In this 
paper, for traditional purposes we use the second approach. For 
RR Lyrae stars the typical value of this difference is about some 
hundredths of a magnitude (ao <^ 1). It is evident that this con- 
stant differs from the zero point of the light curve Ao given in the 
apparent magnitude scale. 



3.1.1 Sinusoidal amplitude modulation 

In the simplest case the modulation is sinusoidal: 

(7m(t) = am sin(27r/mt + v^m). (21) 

Sample light curves obtained with this assumption from Eq. l |20[ ( 
are shown in Fig. [5] Introducing the modulation depth as h — 
a-m/U* the parameters were chosen as ao ^ L^c and am ^ ao, 
resulting in modulations symmetrical to an averaged value viz. a 
horizontal line (left panels). Right panels show cases with higher 
modulation depths (am > ao), where this symmetry is broken. 
A common feature of these light curves is that the maxima and 
minima of the envelope curves coincide in time. Furthermore, the 
average brightness of all light curves vary with /m. It can be seen 
directly from Eq. l |20[ (: the mm(i)ao term is responsible for this be- 
haviour. That is, the found mean brightness {V) variations during 
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Figure 3. Artificial light curves with a sinusoidal AM computed with the 
formula Eq. \20] . Left panels show symmetrical modulation (Om ^ ao; 
ao = 0.2), the right ones are asymmetrical (am > ao; = 0.005). 
The modulation depth h is increasing from the top to bottom as /i = 
0.1, 0.2, 0.4; /m = 0.05 d'^ and (/Jm = 270 deg are fixed. 



Figure 5. Fourier amplitude spectrum of the artificial sinusoidal AM light 
curve in bottom right panel of Fig.|3]after the data are prewhitened with the 
main frequency and its harmonics. Inserts are zooms around the positions 
of the main frequency /o = 2 d~^ (top), and the modulation frequency 
/m = 0.05 d^^ (bottom), respectively. 
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Figure 4. Bottom panel: Artificial light curve with a sinusoidal AM com- 
puted with the formula \20\ . The modulation depth is h = 1.2. Other pa- 
rameters are ao = 0.01, /m = 0.05 d~^ and i^m = 270 deg. Top panels: 
Two-day zooms around a maximum (top left) and a minimum (top right) of 
the modulation cycle. 



the Blazhko cycle ( [Jurcsik et al. |2005) is a natural consequence of 
the AM. 

There is a fascinating case, when the modulation is very strong 
i.e. when the modulation depth is /i > 1. Beside the strong light 
curve changes (Fig. |4]l in some Blazhko phases the shape of the 
light curve looks very unfamiliar (see top right panel in Fig.|4]l. The 
relevance of this mathematical case is corroborated by the Kepler 
observation of V445 Lyr that shows similar characteristics (fig 2. in 
IBenko et al. pOTOt - 

Using some trigonometrical relations, Eq. \2Q) with can 
be converted to a handy sinusoidal decomposition form from where 
the Fourier spectrum is easily seen: 



AM (^) = «o + hai) sin {2-Kfmt 



+ </5m) + 



a J sin {2njfot + ipj) + 

n 

■ aj {sin [27r (j/o - fm)t + ipj] + 
i=i 



[2^07o + /m)t + ^+]}, (22) 



where ipJ = ipj — ipm + tt/2, vp+ = (pj + (fm — vr/2. The 
Fourier spectrum of such an AM signal is familiar for RR Lyrae 
stars' experts (see also Fig. [5]l. It consists of the spectrum of the 
non-modulated star (third term) as in Fig. |2] and two equidistant 
side peaks around each harmonic (last term). The amplitudes of the 
pairs of side peaks are always equal: A{jfo ± /m) ~ ajh/2. The 
second term in Eq. \22) , that causes the average brightness varia- 
tion, produces the frequency /m itself in the spectrum. The Blazhko 
modulation frequency is always found in observed data sets ex- 



tended enough (see e.g. |Kovacs 


1995 Nagy |1998||Jurcsik et al. 


20051 20081 IChadid et al. |2010| 


Poretti et al. |2010||Kolenberg et 


al. |2011|l. 



It is a long-standing question whether there is any Blazhko 
phase where the modulated light curve is identical to the unmod- 
ulated one (see |Jurcsik, Benko & Szeidl ||2002] and further refer- 
ences therein). In this simplest case the answer is easy. It happens 
if the second and fourth terms in \22) disappear simultaneously, 
namely, in the zero points of the modulated sinusoidal function at 
t — {kn — (prn) / {2tt fm) , k is an arbitrary integer. 

The number of used parameters for solving such a light curve 
(Fig. [5} in the traditional way (according to Eq. [T9j is 6n + 5, 
where n denote the number of detected harmonics including the 
main frequency. The necessary number of parameters in our han- 
dling is 2n + 5. The modulation is described by 3 parameters (/m, 
dm, fm) as opposed to the traditional framework where this num- 
ber is 4n + 3. 
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3.1.2 Non-sinusoidal AM 



As a next step, we assume the modulation function m^{t) to be 
an arbitrary periodic signal represented by a Fourier sum with a 
constant frequency /m- Substituting it into Eq. pOl we get 



nilM(i) 



p=i 



A ■ 

ap sm 



J- , A 



c{t), (23) 



where constants are defined by a* = 1 + (a™/[/c ), and = 
/U* . From now on, the upper index A denotes the AM parame- 
ters. Some typical light curves are shown is Fig.|6] It is evident that 
their envelope curves are non-sinusoidal and their shapes depend 
on the actual values of and ip^. The maxima and minima of 
these envelope curves occur, however, at the same Blazhko phase 
as in the previous sinusoidal cases. Rewriting of l |23| l similarly to 
Eq. \22) but in a more compact form yields 



p=0 j=0 



[2^ (j,fo±p/m)i + ^fp] 



(24) 



where the two sinusoidal terms appearing analogously to Eq. \22) 
are formally unified into one formula and denoted by ± signs; 

= "^J +'Pp - 7^/2; f]p = "y^j - "^p + 1-/2. The arbitrary 
constants are chosen to be ipo :— ipo :— tv/2. 

By investigating the Fourier amplitude of the side peaks we 
found that ^(j/o ±p/m)/^(j/o) ~ a^- (i) The amplitude ratio of 
the side peaks of a given order and the central peak is constant, (ii) 
The commonly used amplitude ratio A{j fo±pfm) /^(,/o±P,/m) ~ 
aj/ai vs. harmonic order is the same as the amplitude ratio of the 
main frequency A{jfo)/A{fo) ~ aj/ai vs. harmonic order (iii) 
Since the same coefficient belongs to both side peaks (at ±p), 
the amplitudes of left-hand- side and right-hand-side peaks are the 
same. According to this, the generated Fourier spectrum (Fig|7| 
now shows symmetrical multiplet structure of peaks around the 
main frequency and its harmonics (jfo ± pfm). Each multiplet 
structure is the same at each harmonic order, that is the number 
of the side peaks, their frequency differences and amplitude ratios 
to their central peaks are the same. It is important to note that the 
number of side peaks (in one side) is equal to p. In addition, the har- 
monic components of the modulation frequency p/m also appear. 
(This can be obtained from Eq. ([24j if j — 0.) 

Such a phenomenon was undetected in the observed data of 
Blazhko stars until recently. |Hurta et al. [(2008 1 found equidistant 
quintuplets in the spectrum of RV UMa for the first time. Besides 
triplets and quintuplets, sextuplet structures were also found by|ju 



rcsik et aT"] ( |2008^ in the spectrum of MW Lyr, while [Chadid et 



al. 1 2010b detected even 8th order (sepdecaplet) multiplet frequen- 



cies in the spectrum of CoRoT data of V1127Aql. According to 
Sec. |2.2[ the angle modulations cause infinite numbers of side peaks 
around each harmonic, therefore, the origin of the observed multi- 
plets as a non-sinusoidal amplitude modulation is certain for those 
cited cases (e.g. MWLyr and V1127 Aql), where the harmonics of 
the modulation frequency are also detected. 

In searching for a Blazhko phase where the modulated and 
carrier waves are identical we concluded that the modulation 
terms can only be entirely disappearing from the formula l |24[ ( if 
fflo^ = 1 (a™ = 0) is true, otherwise no such Blazhko phase ex- 
ists. This necessary condition is complemented by an additional 
one: the modulation virtually disappears in the moments when 
J2l^i "jCtp sin[27r(j/o ± p/m)i + ^pfp] = 0. The sum has either 
zero or infinite numbers of zero points depending on the values of 
the parameters a^, ipp. That is, generally there are no such phases 
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Figure 6. Synthetic light curves of non-sinusoidal AM signals computed 
by the formula \23\ . A two-term sum of modulation signal was assumed: 
= 0.01, = 0.2 mag; ip^ = 270 deg are fixed and the phase of 
the second modulation term varies from top to bottom as: = HO. 140, 
220, 270, 360, respectively. 
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Figure 7. Fourier amplitude spectrum of the synthetic non-sinusoidal 
AM light curves in Fig. |6] after the data were prewhitened with the main 
frequency and its harmonics. Inserts are zooms around the positions of 
the main frequency /o = 2 d~^ (top), and the modulation frequency 
/m = 0.05 d^^ (bottom), respectively. 
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where a non-sinusoidal AM light curve and its carrier wave are 
identical. 

The necessary number of the parameters for a light curve fit of 
( [19} and (|23| is {2q + l)2n + 2g + 3 and 2n + 2g + 3, respectively. 
Here n denotes the total number of used harmonics including the 
main frequency and q is the order of side peak structures as above, 
(i.e. q = ^ means triplets, q = 2 quintuplets, etc). In the traditional 
description each additional side peak order increased the number 
of terms by 4n + 2 as opposed to our method, where this increment 
is only 2. 

3.1.3 Parallel AM modulation 



Multiperiodic modulation was suspected in XZ Cyg fLaCluyze et 
al. 2004 ), UZ UMa ( |S6dor et al. |2006ll, SUCol (Szcz ygiel & Fab^ 
rycky ||2007| , and LS Her jWils, Kleidis & Broens ||2008| l. The 
Blazhko RR Lyrae stars of the MACHO and OGLE surveys JaT 



|cock et al. |2000[ [Moskalik & Poretti |2003) that have unequally 
spaced triplet structures in their Fourier spectra are possibly also 
multiperiodically modulated variables. CZ Lac ( |S6dor et al. |201 1} 
is the first Blazhko star with multiperiodic modulation where both 
modulation periods are identified. Not only modulation side peaks 
but linear combinations of the modulation frequencies also appear. 
Signs of multiple modulation were discovered in Kepler data of 
V445 Lyr (Benko et al. 2010). There are numerous possibilities for 
creating a multiply modulated light curve. Let us review some of 
them. 

The most simple case is a natural generalization of Eq. p.3\ 
when the modulation signal is assumed to be a sum of signals with 
different where r — 1,2,... signs constant frequencies. Let 
signals be independent, i.e., the modulation signal consists of lin- 
early superimposed waves. In this case, Eq. \23) reads as: 



»TlAM(i) = 



r—l p—l 



t + ippr 



where ao = 1 + X^^^i o-oi-fU*, and cipr = a^^/U^ 



(25) 

. This formula 



is demonstrated in Fig. [8] In this figure, only two modulation waves 
are taken into account and the only varied parameter is the fre- 
quency of the second modulation is When the two modulation 
frequencies are comparable (/^ = 0.1 d^^ and f^^ = 0.09 d^^ in 
panel a the envelope shape of the light curve shows the well-known 
beating phenomenon. Here the beating period is 200 days, albeit the 
modulation periods are close to the shortest known ones. It is easy 
to understand that the observations taken on a moderate time span 
often detect only the gradual increase or decrease of the amplitude 
of the Blazhko cycles. In panel b of Fig. [i]/?^ was set to 0.075 d"\ 
that is the ratio of the modulation frequencies is 4:3, similarly to 
the case of CZLac during its second observed season in Sodor et 
|al. | ( [201l) . The amplitude changes of the consecutive Blazhko cy- 
cles need well-covered long-term time series observations, other- 
wise the interpretation becomes difficult. Panel c in Fig. [8] shows 
the case where the frequency of the second modulation is half of 
the first one (/^ = 0.05 d~^). These specially selected values 
cause alternating higher and lower Blazhko cycles. The exact 2:1 
ratio between the two modulation frequencies leads to the same re- 
sult as a two-term non-sinusoidal modulation in Eq. \23\ (see also 
top panel in Fig.[6](. The bottom panel in Fig.[8]shows a case where 
the second modulation has a much longer period than the primary 
Blazhko cycle. In a first inspection the top and bottom panels are 
very similar apart from a phase shift. 

To reveal the real situation we need to compare their Fourier 
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Figure 8. Artificial light curves calculated with two independent sinusoidal 
AM modulations according to the formula of (25). The fixed parameters 
wereao = 0.01, an = 0.5, ai2 = 0.2 mag, O-ld^^i^n = 270, 
1^12 = 120 deg, where changes from top to bottom as 0.09, 0.075, 0.05 
and 0.01 d~^, respectively. 



spectra starting with 



= 1 p=0 j=0 



J- I -± 

t + ipjp, 



Where the constants are chosen similarly to j24|: 'Pjpr 



(26) 

^ ^ - 

ifipr + 7r/2; (fi^p^ — (fij + 0pr — 7r/2; ipo — (for = 7r/2. It is easy 
to see that the Fourier spectrum of expression \26\ contains s sets of 
side peaks shown in Fig. |7] The qualitative structure of these sets 
is the same. It consists of the carrier's spectrum (j/o), the peaks 
of the different modulation frequencies and their harmonics (pf^), 
and the side peaks around the main frequency and its harmonics: 
jfo ± p/m, where p — 1,2, ... , q,., and r = 1, 2, . . . , s. Due to 
the independence of the modulation waves no further side peaks 
appear. 



3. 1.4 Modulated modulation - the AM cascade 

It is hard to imagine, however, that in a real stars's case, the dif- 
ferent modulating waves are independently superimposed without 
any interactions. Let us investigate the possibility of the modulated 
modulation: the cascade. In other words, the modulation signal is 
composed of recursively modulated waves as 
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Figure 9. Fourier amplitude specfra of tlie artificial cascade AM light curves 
with parameters of panel b in Fig. [8]- (black) dotted line - and its cascade 
equivalent - (red) continuous line. The spectra show the interval around the 
main pulsation frequency after the data are prewhitened with it. The (blue) 
dashed line shows the spectrum of cascade case after the side peaks /oit/,Jj, 
/o it are also removed (for better visibility the top and bottom spectra 
are shifted with +0.01 and —0.01 mag, respectively). 



„(i) 



{t) ~c'{t), mii'(f) =m;,(t), 



(t) m^ll,{t) = [l + m'^>{t)]c^'>{t), 



As) 



{t) = [l + m[:\t)]c''>{t). ill) 



(s)i 



■m*AM{t) = n 



ao, 



+ ^ d.pr sin \2-KvflJ' + 



tpp,. 



c(t), 



(28) 

where ao,. = 1 + aol./U*,., d,pr = a^r/U'r- U^r denotes the am- 
plitude of the rth carrier wave c'*^' (t). On the basis of a visual in- 
spection, there are imperceptible differences among the light curves 
produced by this expression \2S\ and those that can be seen in 
Fig. [S] The Fourier spectrum, however, contains additional peaks 
at the linear combinations of fo and as it is shown in FigM To 
understand this spectrum we generate Eq. ([28j in a form similar to 
Eq. (I26). 



S Qr n ^ 



I 



\ 



r=l p=0 j=0 I T.' = \ 

sin (^jfo ± j3,/m ]t + ip 



I] I] n ''Or 

fc=2rG7JfcP = l r'STCg 



=1 



fc=2 reTCt 



r'eTZ): 



iP), (29) 



Vjpr = V^J - "fpr + ^/2; (P'jpr = ] + <^pr " 7r/2; (por = 7r/2. 

An index set TZa is defined so that it contains all r indices from r — 
1,2, ... ,s. Index sets TZt means all subsets of 7?.s what contains k 
elements. Therefore, the total number of TZk sets is (^) . The sums 



over r £ TZk mean a sum over the all possible combinations of the 
k number of different indices r. Similarly r' always runs over the 
complement of a set TZk of the actual TZk. The functions Sk build 
up sums of sinusoidal functions (see for the definition and further 
details in Appendix|Bj of the linear combinations of k angles. Here 
the components of the a and /3 vectors are ar — 2TTpf^t + ippr, 
r e TZk, Pk+i = 2njfot + If J and Pr = a^. 

For better comparison we present the formula in the form 
of Eq. \29\ instead of the most possible compact one. Although 
Eq. l|29| seems to be complicated, the meaning of each term is 
simple: the first term can be directly compared to the linearly su- 
perimposed case Eq. \26\ . It produces all the peaks in the Fourier 
spectrum appearing in that case: the main pulsation frequency and 
its harmonics jfo, modulation frequencies and its harmonics p/m, 
the side frequencies around the main frequency and its harmonics 
(jfo ± P,fm)- The second sum in Eq. l |29| belongs to the peaks of 
all the possible linear combinations of p/m, whilst the last term 
is responsible for the many linear combination frequencies around 
the main frequency and its harmonics ( Fig. [9}. The la tter two types 
of combination terms were detected by |S6dor et al. |p01I| > in the 
Fourier spectrum of CZ Lac, the only well-studied multiply modu- 
lated RR Lyrae star. 

Long-term secondary changes in Blazhko cycles can be ex- 
plained by the variable strength of the modulation. To formulate 
this assumption we arrived at 

mlM W = {1 + [1 + mUt)] m[Ut)}c*{t) (30) 

The formula l |30^ can be considered as a special case of (|28j when 
s — 2 and Soi — 0. 



3.2 Blazhko stars with FM 

We remind the reader of the possible absence of real Blazhko stars 
with pure AM that was mentioned in the introductory paragraphs 
of Sec. |3.1| The only difference between pure AM and FM cases is 
that RR Lyrae stars showing pure PM/FM are much more rarely re- 
ported than pure AM ones, but there are some examples (e.g. |Kurtz| 
let al. |2000||Derekas et al. |2004| l. 

How can the formalism discussed in Sec. |2.2| be applied to 
RR Lyrae stars? Let us assume the same carrier wave as in the case 
of AM, but here the instantaneous frequency f{t) is denoted as /o + 
7711*1 (t), and m^it) is an arbitrary (bounded) modulation signal. 

n 

mimit) = ao + X] aj sin{27rj [fo + ml,{t)] t + ipj}. (31) 

j=i 

Expression jSTJ describes a general frequency modulated RR Lyrae 
light curve. 



3.2.1 The sinusoidal FM 

When the modulating function is sinusoidal and expressed in the 
same form as |2TJ, Eq. |3TJ becomes 

»T^FM(i) = ao + 

n 

aj sin |^27rj7ot + jcf sin (2tv f^t + ip^^ + i^j j , (32) 

wherea^ = am/.fui,^^ ~ <^m + 7r/ 2 and the upper index F marks 
the parameters of FM. The amplitude of this signal is determined by 
the Fourier amplitudes aj of the carrier signal, hence no amplitude 
changes are present. In the bottom panel of Fig.fTOla simulated light 
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Figure 10. Bottom: Artificial FM light curve produced with sinusoidal 
modulation according to Eq. (321. Fourier parameters of the "carrier light 
curve" are the same as before, a'' = 0.279; <fi^ = 0. Boxes show the lo- 
cations of the top panels. Top panels: One-day zooms from two different 
phases of the modulation cycle. The non-modulated "carrier" light curve 
is shown by a (red) continuous lines while (blue) dotted lines denotes FM 
signal. The periodic phase shift (viz. PM) caused by FM can be clearly 
identified. 
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Figure 11. Bottom: Fouiier amplitude spectrum of the artificial sinusoidal 
FM light curve in Fig. |10| after the data are prewhitened with the main fre- 
quency and its harmonics. Top panels are zooms around the positions of the 
main frequency /o = 2 d~^ (top left), and its 7th harmonics 8/o = 16 d~^ 
(top right), respectively. The modulation frequency fm is missing from the 
spectrum (insert). 



curve is shown. It is evident that there is no amplitude change. Two- 
day zooms from two different phases of the modulation cycle are 
shown in the top panels. The periodic phase shift caused by FM can 
be identified well: in the left-hand-side panel the non-modulated 
light curve is to the left from the FM light curve whilst in the right- 
hand side panel the situation is opposite. 

Using Chowning's relation \12\ we get from ^32\ : 



n oo 

ajJk (j"-^^ sin ^2tv {jfo + kf^) t + kip^ 



j=i fe=- 



(33) 



This equation shows the main characteristics of the Fourier spec- 
trum (Fig. [TTJ. It consists of peaks at /o and at its harmonics 
jfo and each of them is surrounded by side peaks at jfo ± fc/m 
with symmetrical amplitudes at the two sides. This symmetry of 
the amplitudes can be seen from the expression of amplitude ra- 
tio A{jfo ± kfm)/A{jfo) ~ \J±k{ja^)\, and it is known that 
J-k{z) ~ {—l)^Jk{z). It is worth to compare the AM spectra 
in Figs. |5] and|7] to this FM spectrum. The Fourier amplitude of 
the side peaks are proportional to the Bessel function, and an im- 
mediate consequence can be seen in the figure: the amplitude of 
the triplet peaks are higher at 3/o than at 2/o. (Although it is not 
shown in the figure, the higher order j > 5 harmonics have also 
smaller amplitudes than their side frequencies' ones.) Since the ar- 
gument of the Bessel functions depends on the order of harmonics 
j, higher order harmonics "feel" larger modulation index, which 
results in more side peaks around the higher order harmonics (see 
inserts in Fig. |1 l[l. This effect was found for V1127Aql from its 
CoRoT data by Chadid et al. | ( [20I0| |. A further remarkable flavour 
of this Fourier spectrum is, that it does not include /m as opposed 
to any of the AM spectra (insert in Fig.[TTJ. 

Let us return to the question whether there is any phase in the 
modulation cycle where the light curve is identical to the monope- 
riodic light curve (carrier signal). As it was shown in the case of 
sinusoidal AM modulation, there are some such possible phases. 
Looking at the formula (|32|l it can be realised that the modulation 



disappears at the moments of time if t = {In — ip )/(27r/ni), where 
I is an arbitrary integer 

Estimating the number of the necessary parameters for a light 
curve fit the traditional description Eq. \19\ needs ~ 2n + 3 + 
4 X]j=i [iiit(ja^) + 1], where "int" means the integer function, and 
n is the number of all harmonics including the main frequency as 
well. At the same time, Eq. l |32| ( requires only 2n+5 parameters, no 
more than in the sinusoidal AM case. For a typical case plotted in 
Fig. 10 (n — 10 and = 0.27), the difference is 143 parameters 



versus 25. 



3.2.2 The case of non-sinusoidal FM 



Assuming an arbitrary periodic modulation with a fixed frequency 
we substitute a Fourier sum representing this modulation signal into 
Eq. ^3l\ and get 



n^FM(i) ~ ao + ^ aj sin 



27rj/ot 



j ^ Op sin [2Tvpfmt + iflj + < 



, (34) 



where the constant terms are contracted as (f>j ~ jag + ifj. As in 
the previous sinusoidal case the equation can be rewritten as 



»T^PM(i) = "0 + 



j — 1 ki ,k2 , . ,kp — — oc 

1 \ 1 



n ■^''P (■^^ 
.P=l 



271" ( jfo + ^ kppfm U + X] kpifil + , 
p—1 / p—1 



(35) 
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Figure 12. Bottom: Fourier amplitude spectrum of an artificial non- 
sinusoidal FM light curve calculated by the formula of |34j after the data 
are prewhitened with the main frequency and its harmonics. Parameters of 
the generated light curve were the same as for the light curve in Fig.Qo] and 
p = 2, = —0.1 mag, (/Jj = 7r/4. Top panels are zooms around the 
positions of the main frequency /o = 2 d~^ (top left), and its 7th harmonic 
8/o = 16 d~^ (top right), respectively. 



In one sense, this formula is a generalisation of \\2)) to the case 
for a non-sinusoidal carrier wave, in the other sense however, the 
modulation frequencies are chosen specially as := p/m- 

Comparing equation \35\ with |33j it can be realised that the 
structure of both Fourier amplitude spectra should be similar (cf. 
also Fig.[TT]and Fig.|12[l, although there are significant differences, 
as well. First of all, besides the same values of common Fourier pa- 
rameters of a sinusoidal and a non-sinusoidal case, the detectable 
side peaks are more numerous than for the non-sinusoidal case. The 
reason is simple: the higher order terms in the modulation signals' 
sum increase the "effective modulation index". The most notewor- 
thy difference is the disappearance of the symmetry between the 
amplitude of the side peaks from the lower and higher frequency 
parts. 

To understand this let us investigate the simplest non- 
sinusoidal case, if g = 2 and concentrate only on the side peaks 
around the main pulsation frequency (j = 1). Then the above ex- 
pression \2i5\ is simplified to 



oo oo 



— — oo ^2 — '"'^ 

|27r [/o + (fci + 2k2) U] t + fclV9i + fc2<y32 + <?!>i} ■ (36) 



For calculating the amplitude of the triplet's side peaks A{fo± fm) 
we have to sort out the corresponding terms from the above infinite 
sum such as fci = 1 — 2^2 and fci — — (2fc2 + 1), for the right- 
hand-side and left-hand-side peaks, respectively (k2 is an arbitrary 
integer). It can be seen, that both sums include the same elements, 
because J_3(af)Ji(a2) = J3(af ) J_i(a2 ), J-sloi) ^2(02) = 



for each pair and the relative phase differ- 



ences have the same values with an opposite sign. The only dif- 
fering terms contain Ji(ai ) Jo(a2 ) in the sum of A{fo + fm) and 
J_i(af ) Jo(a|) in the sum of A{fo — fm)- These terms are re- 
sponsible for the asymmetry of the side peaks. Introducing power 
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Figure 13. Amplitude ratios of the harmonic components of the main pul- 
sation (A{j fo) /A{fo)) plotted on decimal logarithmic scale compared to 
the amplitude ratios of the modulation components: A{j fo+pfm) /A{fo + 
Pfm)- Top: sinusoidal FM case; p = 0, 1, 2 and 3 (black) asterisks, (red) 
circles, (blue) triangles and (green) squares, respectively. Bottom: non- 
sinusoidal FM; shape of the symbols denotes p as in the top panel, but 
filled symbols mark positive ps (higher frequency side peaks) while the 
open symbols show negative ps (lower frequency side peaks). 



difference of the side peaks as in Sec.|2.3|we get 



(37) 



Here Ai and $1 indicate the amplitude and phase of a sinusoidal 
oscillation obtained by summing all the terms in ([36j except the 
different ones. The asymmetry of the higher order side peaks (| fci + 
2fc2j > 1) can be verified in a similar manner. 

This asymmetry has a further consequence. The functions of 
amplitude ratio vs. harmonic orders belonging to a given pair of 
side peaks are diverge from each other (Fig. |13^ . This behaviour 
is well-known from the similar diagrams of observed Blazhko 
RRLyrae stars (|Jurcsik et al. |2009b; Cha did et al. |'2 010; Kolen^ 



|berg et al. |20TT 



. It can also be seen that the actual character of the 
asymmetry can change with the harmonic order j or even within a 
given order with the different p. For example, in Fig.[T3]if p = 1 
(triplets), the right-hand-side peaks are always higher than the left- 
hand-side ones and the difference between the pairs are increas- 
ing with harmonic orders. Meanwhile, if p = 3 (septuplets) the 
situation is the opposite. In the case of p = 2 (quintuplets) the 
lower frequency peaks have higher amplitude around the lower or- 
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Figure 14. O-C diagram of the maxima of FM light curves with sinusoidal 
- (blue) filled circles - and non-sinusoidal - (red) empty circles - modu- 
lations, respectively. The input light curves are generated from the formu- 
lae ofp2j and J34| with the same parameters as the light curve shown in 
Fig. 1 1 0| (sinusoidal case) and Fourier plot in Fig. |12| (non-sinusoidal case). 
(For better visibility the non-sinusoidal curve is shifted by —0.05.) 



Figure 15. Fourier amplitude spectra of the artificial light curves containing 
two parallel FM modulations in formula J38}. The figure shows a zoom 
around the main pulsation frequency after the data are prewhitened with it 
and the triplet components /o it f^-^, /o it f"^-^. The highest peaks are at the 
quintuplet frequencies (/o it 2/j5j, /o it 2/,^) and at the hnear combination 
frequencies. (The used parameters are: = 0.1 d~^, = 0.01 d~^, 
a\-^ = 0.5 mag, = 0-2 mag.) 



der (j < 5) harmonics, for higher harmonics (j > 7) the amplitude 
ratios of the pairs of side peaks reversed. 

As it was discussed in the introductory Section |2.3| simulta- 
neous and sinusoidal amplitude and phase modulations result in 
an asymmetrical spectrum, therefore, the asymmetry of the ampli- 
tude spectrum alone is not a good criterion for detecting a non- 
sinusoidal FM. The classical O— C diagram is an ideal tool for this 
purpose (see e.g. |Sterken |2005p . Fig. |14| illustrates the O— C dia- 
grams of the maxima for two artificial FM light curves: a sinusoidal 
and a simple non-sinusoidal one. 

At the end of this section we compare the necessary param- 
eters of a potential fit based on the classical description ljT9j and 
the present ^3l\ one. In the latter case this value is 2n + 2q + 3, 
where n and q are defined in ^3l\ . The expression is the same as 
in the case of a non-sinusoidal AM. The traditional formula needs 



^ 2n + 3 + (4E; 
case showed in Fig. 
values are 27 and 163, respectively. 



10 



12 



(n 



+ 1]) parameters. For the 
= 0.27, = 0.1) these 



3.2.3 Parallel FM 

We continue the discussion as in the case of AM. The next step 
is the multiply modulated FM with independently superimposed 
modulation signals (parallel modulation). As has already been 
noted, the chance of such a scenario is very low for stars, but this 
case shows a new phenomenon, which is why it is worth to have a 
look at it. 

r 

rhFyiit) = do + ^ cij sin < 2njfot + 



^ ^ Qpr sin hivpf^t + (fil. 



r — l p—l 



+ > 



(38) 



Here < 



j&Q^ + ipj. The formula 1 38 i can easily be transformed 



with the help of Eqs. (jTSj and l |35| l to 



= 1 feii,fei2,....fe,gs = 



r — l p—l 



r—l p—l 



(39) 

There is a fundamental difference between the construction of 
Fourier spectra of parallel AM and FM signals. While the AM spec- 
tra build up from a simple sum of the component spectra belonging 
to a given modulation frequency FM spectra contain all possi- 
ble linear combinations of and the harmonics of the main pul- 
sation frequency j/o. This is illustrated in Fig.[T5] In practice, this 
effect complicates distinguishing Fourier spectra from the cascade 
AM and parallel FM. 

3.2.4 The FM cascade 

Although a parallel FM modulation results in a more complex 
Fourier spectrum than either a parallel or even a cascade AM, our 
former statement is still true. There is a very low chance for in- 
dependently superimposed modulation signals in real stars. Let us 
turn to the FM cascade (viz. the modulated modulation) case now! 



where 

Cfm it) 



= 0,0 



n 

+ ^ a, sin \2-Kjfot + jCym (t) + 4>j] , (40) 



"t-fm(*) = ^ «pi sin [27rp/it + pm'pl^{t) + i^^ij 

3=1 
Is 

mpliit) = ^ fflps sin h-Kpf^t + <^pj , (41) 
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and (f)j = jSgi + (fj. As we can see, the function CFM{t) consists 
of a modulation cascade witli s elements, where all elemental mod- 
ulation functions nip^ (t) are represented by finite Fourier sums. 
That is, they are assumed to be independent periodic signals with 
the frequencies Since an FM modulation can be reproduced by 
infinite series of sinusoidal functions (see Chowning relation), it is 
not a surprise that the sinusoidal decomposition of the expression 
IpOl is very similar to the parallel case \39\. Namely 



J — 1 fcll,fcl2,...,fcg. 



-F 



,p=l 



r=2 p=l 



\ r — l p—l 



(42) 



The frequency content is exactly the same as in the parallel case, 
only the values of amplitudes and phases are different. 



3.3 The case of PM 

Here we discuss the phase modulation. As we stated in the Sec. |2] 
there is no chance to distinguish between FM and PM phenomena 
on the basis of their measured signals (inverse problem), if the mod- 
ulation function m'^{t) in Eq. \3i\ is allowed to be arbitrary. At the 
same time, if the basic physical parameters such as effective tem- 
perature, radius and log g are changing during the Blazhko cycle as 
was found recently ( jSodor, Jurcsik & Szeidl |2009| | Jurcsik et al.| 
[ 2009a b ) the cyclic variation of the fundamental pulsation period 
(vis. frequency) that results in FM would be a plausible explana- 
tion for observed effects. There is an additional possible argument 
against the existence of PM in RRLyrae stars. 

If we assume that the modulating function m*^ contains no ex- 
plicit time variation - as in the usual definition for PM modulation 
in electronics - Eq. \3l\ reads as 

n 

'TipM W = ao + ^ aj sin [2TTjfot + m*^{t) + ipj]. (43) 

When this formula is expressed as Eq. or Eq. \35) according 
to a sinusoidal or an arbitrary periodic modulating function, respec- 
tively, the arguments of Bessel functions are independent from the 
harmonic order j as opposed to the case of FM. It causes a sys- 
tematic difference between Fourier spectra of FM and PM. While 
the number of detectable side peaks in FM increases with the order 
of harmonics, for PM the number of side peaks is the same for all 
harmonics. 

There are two Blazhko RRLyrae stars that show both strong 
phase variations and their data are precise enough, these are the 
CoRoT targets V1127Aql and CoRoT 105288363 ( jChadid et al.| 
[] [20T0{ [Guggenberger erar][20TT) . The spectrum of V1127Aql 
clearly shows the existence of FM: 3rd order side frequencies 
are detected around the main pulsation frequency while order of 
8th around the 19th harmonic. The Fourier analysis of separate 
Blazhko cycles of CoRoT 105288363 showed that with the increas- 
ing strength of phase variation, the number of detected side peaks 
around higher order harmonics are also increased jGuggenberger| 
|et al. |20TT| . It is an evidence of (changing) FM. 



3.4 Real Blazhko stars with simultaneous AM and FM 

In this section we discuss the general combined case, when both 
types of modulations occur simultaneously. As it was mentioned 
before both AM and FM type modulations were detected for all 
observed Blazhko RRLyrae stars if the observed data sets were 
precise and long enough. This is the situation for ground-based 
( [Jurcsik et al. |2009c[ and references therein) and space-bom ob- 
servations of CoRoT and Kepler ( Chadid et al. |20I0{|Poretti et al. | 



[20T0) [Benko et al. ]20T0i [Kolenberg et al. [2011^ aswell 

Generalising the sinusoidal case of combined modulation 
Eq. discussed in Sec. |2.3| we get 

"icomb (i) = [1 + (t )] mpM it) , (44) 

where m^M(t) is the general modulated FM function defined by 
Eq. \3\) . Since all observed Blazhko stars show AM and FM with 
the same frequency, we have investigated only those cases where 
this assumption is fulfilled. 

3.4. 1 Combined modulations with sinusoidal functions 

The simplest case similarly to the pure AM and FM cases is the 
simultaneous but sinusoidal modulations. 



n^Comb(i) = (1 + hsm2nfmt) ■ 

au + ^ aj sin |^27rj/ot + ja^ sin (2nfmt + (pm) + (/'jj | , 

(45) 



where the notations are the same or directly analogous with the 
previously defined ones: h — am/U^'M and U^m is the amplitude 
of the second term (the FM modulated "carrier wave"). The relative 
phase between AM and FM signals is = ip^ — ipm- 

According to the schema of \16\ expression ( [45 ^ can be refor- 
mulated into 

mcomhii) = ao + aohsm2-Kfn^t+ 



E E 

J=l k——c 



^Jk{ja^) sin [27r (/o + kfm) t + kcpm + <y5j] + 
{ja') sin [2tt (/o + k,U) t + [k ~ 1) + ^pj] + 



h 



Jk+i{ja') sin [27r (/o + fc/^) t + {k + l)<l>^ + ip-] 



(46) 



where p- = ipj ± 7r/2. Based on Eq. (45 ' the Fourier spectrum 
in Fig.[T6]can be interpreted as a sum of tEe combined modulation 
with sinusoidal carrier wave \i6) with an additional term describ- 
ing the modulation frequency itself (insert in bottom panel). Each 
harmonic is surrounded by a multiplet structure of peaks just like 
the main frequency. The number of side peaks increases with the 
harmonic order j similarly for FM (Sec. |3.2.T) . 

The asymmetrical amplitudes of pairs of side frequencies be- 
longing to a given harmonic can be characterised similarly to the 
sinusoidal carrier wave case jl8[ l as 

hi 



Aj, = -4- 



(47) 



Here Aji = A^ijfo + - A^{jfo - l.fm) is the power dif- 
ference of the Ith side peaks at the jth harmonics (/ — 1,2,...). 
Similarly to the course book case discussed in Sec. [2.3[ the asym- 
metry depends on the actual value of h and , (viz. the relative 
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Frequency [d ] 

Figure 16. Bottom: Fourier amplitude spectrum of an artificial com- 
bined (AM & FM) light curve computed from Eq. \A5) after the data are 
prewhitened with the main frequency and its harmonics. Top panels are 
zooms around the positions of the main frequency /o = 2 d~ ^ (top left), 
and its 7th harmonics 8/o = 16 d~^ (top right), respectively. The relative 
phase between AM and FM is set to 0m = 270 deg. 



strengths of AM and FM) and the relative initial phase angle <j)m- 
The most extreme possibility is when one of the side peaks com- 
pletely disappears. The necessary conditions are = ±7r/2 and 
jcf = hi. The asymmetry decreases with the increasing harmonic 
order j (see also top panels in Fig. [16}, because all the Bessel func- 
tions quickly converge to zero with increasing arguments, therefore 
dominate the right hand side of expression |47]l. 

Non-equidistant sampling and large gaps in the observed time 
series can cause significant differences between side-peak ampli- 
tudes (see |Jurcsik et al."||2005| l. Such sampling effects, however, 
can not explain huge differences, such as when side peaks com- 
pletely disappear in one side and the spectra show doublets, though 
numerous examples were found by large surveys as MACHO and 
OGLE l ,Alcock e t al. |2000||2003||Mosk alik & Poretti |20 03l. But 
as illustrated by Fig. [T6] highly asymmetrical side peaks can eas- 
ily be generated by \A5\ . This asymmetry effect can be a possible 
explanation for the observed doublets (RR-i^l stars) and even for 
triplets (RR-!^2 stars). In the latter case the two side frequencies 
can originate from a quintuplet structure (equidistant triplet on one 
side) or from a multifrequency modulation (non-equidistant triplet 
on one side). 

Searching for phases where the modulated and non-modulated 
light curves are identical we conclude that such phases exist only 
if 4'm = (fci — ^2)71"; iki, k2 are integers) and then the moments 
of the coincidences are t = fc2/2/m. Assuming all Blazhko stars 
showing both AM and FM this conclusion supports the finding of 
|Jurcsik, Benk6 & Szeidl | (|2002), who studied light and radial ve- 
locity curves of Blazhko RR Lyrae stars. 

The amplitude ratio vs. harmonic order diagrams show sim- 
ilar shapes and relative positions as were discussed in Sec. |3.2.2| 
with the connection of Fig. [13] Let us look at the maximum bright- 
ness vs. maximum phase diagrams, a classical tool for analysing 
Blazhko RR Lyrae stars. Such diagrams are plotted in Fig. |17|for 
synthetic light curves generated from the formula of \A5\. All 
diagrams have a simple round shape. They reflect the relative 
strength of AM and FM components. In panels A and B the rel- 
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Figure 17. Maximum brightness vs. maximum phase diagrams for some ar- 
tificial light curves with combined sinusoidal modulations. Between panels 
A and B the relative strengths of AM and FM are changed as (A) h = 0.1 
and aF = 0.2 and (B) h = 0.2 and aF = 0.1, respectively. From panel 
B to F the amplitudes are fixed and only the relative phase (/>ni is changed 
as shown at the upper left corner in each panel. The arrows in panels C-F 
indicate the direction of motion. 



ative strengths are opposite 2h = aF and h = 2aF , respec- 
tively. As a consequence, the loop is deformed vertically or hor- 
izontally. When the angle differs from the special values of 
ln/2, (I = 0, 1, 2, 3, 4), the axes of the loops are inclined to the ver- 
tical horizontal position. This angle also determines the direction of 
motion. If < (j>ui < tt it is clockwise, whilst if vr < 4>m < 2t7 
it is anti-clockwise. (These conditions are the same as it was found 
by S zeidl & Jurcsik [2 009 for sinusoidal carrier waves.) It is note- 
worthy, that the same ranges of also determine the character of 
power difference of the side peaks: if the right hand side peaks are 
higher than the left hand side ones then the direction of motion is 
anti-clockwise and vice versa. 



3.4.2 Non-sinusoidal combined modulation 

On the basis of the previous sections it is easy to define the light 
curves which are modulated by general periodic signals simultane- 
ously both in their amplitudes and phases: 



"lComb(i) 



'(t) 



(48) 



where the functions m^M (^) arid rrip-f^ (t) are defined by Eqs. p.3\ 
and J34I, respectively. Since the two modulations are described 
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by different functions (they are represented by different order of 
Fourier sums), no such simple relative phase can be defined as (pm 
for the sinusoidal case in Sec |3AT] Therefore, we obtain for the 
mathematical form of such a generally modulated light curve: 



n^Comb(i) 



ao + 



V sin (2TTp'fmt + ipp'^ 



ao 



27vjfoU 



j ^ flp sin (2TTpfmt + (fl^ + < 



p=i 



(49) 



where the notations are the same as in Eqs. |23j and l |34| l. This 
expression describes all the discussed phenomena of a light curve 
modulated regularly with a single frequency /m- The envelopes 
of these light curves are very similar to the envelopes of non- 
sinusoidal AM light curves shown in Fig.|6] The light curves show 
non-sinusoidal phase variation as well (see also Figs.|10|and|14[l. 

As in the former simpler cases, the Fourier spectrum can also 
be constructed analytically with the help of the sinusoidal decom- 
position of ^49): 



.it) 



A 



ttoip' sin (2-Kp' fuit + + 



n q 



j — l p'— fci,fc2,....A;g — — oo 

sin { 2iT j 



n -^^p^^ 

.p=i 



p=i 



Here the ^p^^,. ~ kp(pp ± (pp + (t>j T 7r/2, the arbitrary 

constant is chosen as ipQ := 7r/2. The qualitative structure of this 
spectrum is simple and well-understandable on the basis of the pre- 
viously discussed cases. The second term is responsible for the 
appearance of the modulation frequency and its higher harmon- 
ics (see insert in Fig. |7|. The next (infinite number of) terms de- 
scribe a spectrum which is similar to the non-sinusoidal FM spec- 
trum (Fig.|12|l but it also shows the AM splitting which is present 
in the sinusoidal combined case. These effects make the calcula- 
tion of the peaks' amplitude complicated. The asymmetry between 
each pair from a multiplet around a given harmonic is determined 
by two factors: one of them is the the non-sinusoidal nature of the 
FM (Sec. |3.2!2t and the other one is the combination of AM and 
FM (Sec.[3AT}. 

The maximum brightness vs. maximum phase diagrams gen- 
erally show complicated shapes. They could form knots, loops and 
other non-trivial features. A collection of such diagrams is plotted 
in Fig. [Ts] The direction of the motion can arbitrarily change by 
tuning the initial phases. 



3.4.3 Combined multifrequency modulations 

A combined modulation with multiperiodic AM or FM or both can 
be handled analogously to the simpler presented cases. We can 



substitute rn^ii) into the general expression i44i as it was de- 
fined by Eq. (parallel AM) or Eq. l|28]( (AM cascade). Writ- 



ing m^M (t ) as Eq. 1 38 1 (parallel FM) or Eq. 1 40 1 (FM cascade) in 
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Figure 18. Some typical maximum brightness vs. maximum phase dia- 
grams for synthetic light curves with combined non-sinusoidal modulations. 
The relative strength, initial phases and number of used harmonics of AM 
and FM has been varied. 



principle is straightforward. In practice, however, calculating co- 
efficients (amplitudes, phases) is more complicated. The resulting 
light curves and Fourier spectra can be interpreted on the basis of 
their constituents. They do not show new features except their max- 
imum brightness vs. maximum phase diagrams which show time- 
dependent and generally non-closed curves as opposed to those in 
Fig[T8] If the ratio of modulating frequencies are commensurable, 
the curve is closed, otherwise it has a non-repetitive behaviour The 
reason is that if the modulation is described by TV independent fre- 
quencies the proper diagram would be 2Af -dimensional and the 
classical one is only a 2-D projection of it. 



4 PRACTICAL APPLICATION - A CASE STUDY 

To demonstrate how our formalism works in practice, we gener- 
ated two artificial light curves with simultaneous non-sinusoidal 
AM (with two harmonics, p' = 2) and FM (with three harmon- 
ics, p = 3). The light curves are 100-days long and sampled by 
5-minutes in the same manner as the all synthetic light curves in 
the paper We added Gaussian noise to the light curves either with 
rms=0.01 mag (model A) or 10"'* mag (model B), respectively. 
The model A is similar to a good quality ground-based observation, 
while the model B simulates a typical space-borne data set. These 
two artificial light curves (top panels in Fig.[T9j were analysed with 
a blind test (i.e. without knowledge about frequencies, amplitudes 
and phases) both by the traditional way and by our method. 



4.1 Classical light curve analysis 

Constructing the mathematical model Eq. \19) of the light curves 
in the traditional analysis is a successive prewhitening process. It 
consists of Fourier spectra building, fitting the data with the param- 
eters of the highest peak(s) in the spectrum by a non-linear algo- 
rithm, and subtracting the fitted function from the data and so on. 
The process continues as long as significant peaks are detected. At 
the end of this analysis the noise of the residual data reaches the 
observational scatter 
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Figure 19. Artificial liglit curves with Gaussian noise and their Fourier spectra. Model A (ct = 10^^) is plotted in the left, model B (cr = 10^*) in the right. 
The residual light curves and their spectra after the traditional fitting process (middle panels) and the present one (bottom panels). (Residuals and their spectra 
are in different scales for the two models.) 



In the case of model A the highest pealcs belong to the main 
frequency /o and its harmonics. In the further prewhitening steps 
the triplets (/o±/m), quintuplets (/o±2/m), septuplets (/o±3/m), 
nonuplets (/o ± 4/in) and the modulation frequency (/m) were 
found and fitted. The significance level was chosen at S/N — 4, 
where signal-to-noise ratio (S/N) is estimated as [Breger et al. 
( |1993^ . To remove all significant peaks from the spectrum five from 
the undecaplet peaks (/o + 5/m) and two aliases have to be fitted 
and subtracted. (We note that the frequency 2/ni was detectable, but 
under the significance level, therefore it was not fitted.) The resid- 
ual light curve of this process and its Fourier spectrum are plotted 
in the middle panels of Fig. [19] The rms of the residual light curve 
is 10"'^, so we got back the input noise value. The number of fit- 
ted frequencies are 103 and used parameters in this successive fit is 
201. If the frequency of the side peaks are also fitted independently 
("let it free approach") this value increases to 286. 

The process works similarly in the case of model B as well. 
Naturally, many more significant peaks are detectable. From the 
highest peaks to the lowest: /o, its harmonics, the side peaks up 
their orders of six (/o ± 6/m), the modulation frequency and its 
harmonic (/m, 2/ni) are significant. Many alias peaks originating 
from the finite data length are also detectable: 46 such frequencies 
were removed up to the significance level of the sixth order of side 
peaks. We stopped the analysis here at S/N « 40, because we 
already found 168 frequencies and used 368 parameters (or 478 if 
we fit each frequency independently). The resulted light curve and 
its Fourier spectrum are shown in the middle panels of Fig. [19] The 
rms of this light curve is 10"'' an order of magnitude higher than 
the input noise parameter is. 



4.2 Light curve analysis in our framework 

When we apply the approach of this work we need to calculate only 
one Fourier spectrum and a single non-linear fit for each light curve. 
The Fourier spectrum and the characteristics of the light curve help 
us to choose the proper fitting formula and to determine the initial 
values of the fit. 

The Fourier spectrum for any of the light curves A and B pro- 
vide us with the necessary parameters (/o, j, ao, aj, ipj) of the 



carrier wave. The amplitude modulation and its non-sinusoidal na- 
ture in both light curves are apparent. Searching for peaks in the 
low frequency range of the Fourier spectra results in good initial 
values for /m, p — 1, 2, a^, and ipp,. While more side peaks are 
detectable around the harmonics of the main frequency than the 
number of harmonics of /m and the side peaks are asymmetrical, 
an FM has to be assumed. To check its non-sinusoidal nature we 
may prepare e.g. an O— C or maximum brightness vs. maximum 
phase diagram. At the end of this preparatory work we can choose 
the fitting formula l |49| >. We note that to determine the correct initial 
values of and (pp depends on the tool used for finding the fre- 
quency variations. In the worst case, they can be estimated by some 
numerical trials with the non-linear fit. Using initial values that are 
good enough the non-linear fit converged fast for models A and B. 
The algorithm reached the noise levels 10~^ and 10"'* within few 
(less than 10) iterations automatically. The residual light curves and 
their Fourier spectra are plotted in the bottom panels of Fig[T9] The 
number of fitted parameters for both models is only 33. Due to the 
finite numerical accuracy the residual spectra always show a struc- 
ture reflecting the original spectra at very low levels. 

In conclusion our method fits the light curves in a single step 
with much less parameters than the traditional one. In addition we 
avoid the time-consuming alias fitting and subtracting processes. 
The difference in the number of used parameters increases with the 
increasing accuracy of the observed data sets. In our example the 
number of regressed parameters is reduced from 6 times (model A) 
to more than 10 times (model B). Our description has advantages in 
the numerical fit of the ground-based observations as well, but its 
advantages are outstanding in the analysis of the space-born time 
series. 



5 DISCUSSION AND SUMMARY 

In this paper we have investigated mathematical representation of 
artificial light curves. These light curves are defined as modulated 
signals where their carrier wave is a monoperiodic RR Lyrae light 
curve defined by its finite Fourier sum. Different types of periodic 
functions are taken into account as modulation functions from the 
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Table 1. A schematic table for classifying different types of single frequency modulations from their light curves and Fourier spectra. The indices I and k are 
integers, where I is a finite number, while k can generally be infinite. 



modulation 


phenomenon 


AM 


FM 


Combined 


sinusoidal 


amplitude variations: 


yes (simple) 


no 


yes (simple) 




phase variations: 


no 


sinusoidal 


sinusoidal 




side peak structure: 


triplets 07o ± /m) 


multiplets (jfo it fc/m) 


multiplets (jfo it fc/m) 




side peak amplitudes: 


symmetrical 


symmetrical 


asymmetrical 




modulation frequency: 


/m 




/m 


non-sinusoidal 


amplitude variations: 


yes (complex) 


no 


yes (complex) 




phase variations: 


no 


non- sinusoidal 


non-sinusoidal 




side peak structure: 


multiplets (jfo it i/m) 


multiplets (jfo ±kfm) 


multiplets (jfo it fc/m) 




side peak amplitude: 


symmetrical 


asymmetrical 


asymmetrical 




modulation frequencies: 






«/m 



simple sinusoidal to multiperiodic and general ones. The conse- 
quences of these modulation functions and modulation types (AM, 
FM, combined) are reviewed and presented. 

We followed a step-by-step analysis from the simplest case 
to the more complicated ones. The results of this process can be 
summarized as follows: 

• (i) Tuning AM by the used modulation function, namely the 
number of harmonic terms and their amplitudes and phases, the 
synthetic light curves reproduce well the observed shape of the 
Blazhko envelope curves. This is always true for the envelopes 
of maxima, however, the envelope of minima are affected by the 
fine structure of the bump and its phase shift along the light curve 
and changing shapes. These variations in the shock are associated 
to the Blazhko modulation such a way that a displacement of the 
shock forming region occurs over the Blazhko cycle. This induces 
an earlier or later occurrence of the bump in the light curves (see 
|Preston et al. |1965] >. If these variations are strong enough, the en- 
velope curve of minima could be different from our simple mod- 
ulated ones, (ii) We showed that AM with extremely high modu- 
lation depth (h > 1) might be an explanation for the strange light 
curve shape of V445 Lyrae observed by Kepler, (iii) When the AM 
function depends on more than one frequency, (in a parallel or cas- 
cade modulation) we can generate envelope curves with different 
observed phenomena such as beating effect, alternating maxima, 
long-term, periodic amplitude changes. 

• Fourier spectra of AM cases can be easily classified (see also 
Table [T}. For a periodic modulation the signal represented by a fi- 
nite Fourier sum the spectrum of the synthetic light curve shows the 
modulation frequency, multiplets around the pulsation frequency 
and its harmonics, as well. The order of multiplets (the number of 
peaks on one side) is the same as the used number of terms for de- 
scribing the modulation function. As a special case, the sinusoidal 
modulation results in triplet structures and the appearance of /m 
only. For a multiperiodic parallel modulation, the spectrum con- 
sists of a sum of the frequencies of each component modulation, 
while in the case of modulated modulation (AM cascade) spectra 
includes additional peaks at the all possible linear combination fre- 
quencies, as well. 

• For all AM cases the Fourier amplitude of the side frequencies 
are proportional to the given harmonics' amplitude, therefore, the 
numbers and the side peak amplitude ratios compared to the central 
frequencies are the same for all orders of harmonics. Pure AM re- 
sults in symmetrical multiplets: side frequencies of the same order 
have the same amplitudes. 

• The variation of the mean brightness through the Blazhko cy- 



cle is a consequence of the AM modulation. The effect appears even 
for the simplest sinusoidal modulation function, and is closely re- 
lated to the appearance of the modulation frequency itself in the 
Fourier spectrum. 

• Pure FM light curves have no amplitude changes but signif- 
icant phase variations. Their Fourier spectra show multiplet struc- 
ture around the pulsation frequency and its harmonics even for the 
simplest sinusoidal modulation. As opposed to AM the detectable 
number of side peaks is increasing with the increasing harmonic or- 
ders. The modulation frequencies (and their harmonics) are absent 
in the spectra. The side peaks' structure can be symmetrical (for si- 
nusoidal modulation) or asymmetrical (most of the non-sinusoidal 
cases). The non-sinusoidal FM can be characterized by the classi- 
cal O— C tool. The multifrequency FM spectra include all possible 
linear combination frequencies in the parallel case. 

• The Fourier amplitude ratio of the side peaks vs. harmonic or- 
ders are determined by the Bessel functions of the first kind (for the 
sinusoidal case) or their product (for the non-sinusoidal case). As a 
consequence, these relationships are generally non-strictly decreas- 
ing with the increasing harmonic orders. Therefore, we can find 
larger amplitudes for the higher-order side peaks than for the lower- 
orders ones for comparing the same order of side peaks around dif- 
ferent harmonics and different order of side peaks around a given 
harmonic. The latter can be true also for the amplitude of the side 
peaks and the central peaks' (the harmonics') amplitude. 

• The simultaneous AM and FM show all the above mentioned 
effects. The side peaks' amplitudes in the Fourier spectra of the 
combined modulations are generally asymmetrical already for sinu- 
soidal AM and FM. The asymmetry can totally remove side peak(s) 
in one side. The strength of asymmetry depends on the initial phase 
difference of the AM and FM modulation functions. This phase dif- 
ference determines both the inclination and direction of motion of 
the loop in the maximum brightness vs. maximum phase diagrain. 

• In the case of the pure sinusoidal AM and FM (and with 
very special parameters for non-sinusoidal and combined cases as 
well), there are well-defined phases where the modulated and non- 
modulated light curves are identical. However, we have pointed out 
that, when a star shows general combined modulation (even a sinu- 
soidal one), its light curve is always different from any monoperi- 
odic (non-modulated) stars. 

• We showed that our modulation description presented here 
needs much (typically 3-10 times) less parameters for fitting a light 
curve than the classical solution, where all detected frequencies 
have their own independent amplitude and phase. The prize to pay 
for the relatively small number of parameters is the more compli- 



Blazhko RR Lyrae light curves as modulated signals 17 



cated fitting formulae. Our approach fits the fight curves only once, 
instead of the traditional successive prewhitening process, and un- 
affected by most aliasing problems. 

It is equally important to list those observed features that do 
not follow this picture. One of these effects is the phase lag, a dif- 
ference between the Blazhko phases of maxima of maxima and 
minima of minima of the modulated light curves. These phase lags 
can be several pulsation cycles long. In our study, however, all en- 
velopes are highly symmetrical: their maxima and minima are prac- 
tically at the same Blazhko phases. If we see the well-sampled ob- 
served light curves (e.g. CoRoT or Kepler data) showing this phase 
lag, we can realise that the effect builds up from the systematic 
motions and changing the shape of the bump caused by the hydro- 
dynamical shock in the RR Lyrae atmospheres. 

An additional unsolved problem is the distribution of the side 
peaks' asymmetry. The large surveys (OGLE, MACHO) found that 
the spectra of Blazhko RR Lyrae stars more frequently contain 
asymmetrical side frequency patterns, where the higher frequency 
side peaks have the higher amplitude than vice versa (e.g. 74 vs. 26 
per-cents for LMC by jAlcock et al. |2000l l. All formulae suggest 
50-50 per-cents of probabilities, if the distribution of the relative 
phases between AM and FM is uniform. The detected asymmetri- 
cal distribution might have deeper physical origin. 

Our description does not explain the additional frequencies 
such as half-integer frequencies belonging to period doubling, tem- 
porary overtone frequencies and further exotic frequencies in the 



Fourier spectra, recently discovered from the space data i 


Gruber- 


bauer et al. |2007 Kolenberg et al. |2010||Szab6 et al. |2010 


Chadid 


et al. |2010||Poretti et al. '2010; Benko et al. |2010f. Their connec- 



tion with the physics of Blazhko effect is poorly understood and 
their modelling is out of the scope of this paper. 

The main result of this analysis is the demonstration that the 
modulation paradigm gives successfully accounts for most of the 
observed properties of Blazhko RR Lyrae light curves. Moreover, 
this work may help to distinguish between those features that have 
deeper physical origin from those ones which appear simply due 
to the modulation. This framework seems to be very flexible and 
can easily be generalised to a larger fraction of variable stars tak- 
ing into account, e.g., multifrequency carrier waves or stochastic 
modulations. 
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APPENDIX A: ABOUT THE EXACT FOURIER 
TRANSFORMS 

There are many definitions of tlie Fourier transform, differing from 
eachi other by their normalizations. In this paper we adhere to the 
following definition: 

/oc 
git)e-'^^f'dt. (Al) 
- oo 

Since we always transform our equations describing the different 
types of modulated signals to linear combinations of sinusoidal 
functions, it is necessary to know the exact Fourier transformation 
of these elements. 

The Fourier transform of the sinusoidal carrier wave Eq. 
itself is 

T [c(i)] = 7vV2^U, {i cos <p, [S{f-fc)- S{f + /c)] + 

sin<^,[5(/-/c) + 5(/ + /e)]}, (A2) 

where i is the imaginary unit, 5 is the Dirac function. As we see, the 
Fourier amplitude spectrum belonging to c{t) has two components: 
one at positive frequency (centred on +/c) and one at negative fre- 
quency (centred on — /c). We concern only the positive frequency, 
since the negative ones have no physical meaning. 

Demonstrating the use of the above formulae, the positive part 
of the Fourier transformation of the sinusoidal AM signal ^ is 
given: 

-F+ lUAM{t)]=T+ [C{t)] + 
n\j^Um [[C0s(v9c — Ifiin) + i sin(l^c — V'm)] S{f + /c — ,/m) — 
[C0s(v3c + fm) + isin(<^c + </'m)] S{f + fc + /m)] • (A3) 

According to this schema, the investigated signals describing the 
more complicated expressions can also be calculated in a straight- 
forward way. 



Each sin/cos term includes a sum of n terms of the Oi angles inside 
their arguments as ai ± 02 ± 03, • ■ • ± a„. One combination from 
these sets contains I positive and n — I negative angles. The total 
number of such a combination is 2"~^. Each row of the matrix 
Q contains the signs of the angles of a possible set from the total 
number of 2"^^. Namely, the elements of the matrix Q are either 
+ 1 or —1, but Qii — 1 i.e. the first columns contain +1 for all Is. 

This paper has been typeset from a TgX/ KTgX file prepared by the 
author 



APPENDIX B: GENERALIZED PRODUCT-TO-SUM 
FORMULAE 

Let us calculate the generalised form of the product-to-sum expres- 
sion from the well-known identity of trigonometry: 

sin(ai) sin(a2) — ^ [cos(ai — 02) — cos(qi + 02)]- (Bl) 

Assuming n factor in the left-hand- side it can be formally written 
as nr=i sinOi. Introducing the vector-scalar function Sn{cx) : — 
nr=i sinoi, a — (ai, Q2, . . • , a„)'^, where T indicates transpo- 
sition, and applying the formula |bTJ n times recursively we arrive 
to 



S„{a) ~ < 



,(l-n) 



2("-l) 



J2 (-1) cos Qi^ct 



2("-l) 

2(1-") J2 (-1) 



AT' 



J2 Qua, 

i = l 



where 



if n even, 



if n odd, 



(B2) 



N : 



3n 



j^Qu, and N' = ^^^-j2Qi: 



